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Abstract 



An inclusive search for supersymmetric processes that produce final states with jets 
and missing transverse energy is performed in pp collisions at a centre-of-mass en- 
ergy of 8 TeV. The data sample corresponds to an integrated luminosity of 11.7fb _1 
collected by the CMS experiment at the LHC In this search, a dimensionless kine- 
matic variable, otj, is used to discriminate between events with genuine and misre- 
constructed missing transverse energy. The search is based on an examination of the 
number of reconstructed jets per event, the scalar sum of transverse energies of these 
jets, and the number of these jets identified as originating from bottom quarks. No 
significant excess of events over the standard model expectation is found. Exclusion 
limits are set in the parameter space of simplified models, with a special emphasis 
on both compressed-spectrum scenarios and direct or gluino-induced production of 
third-generation squarks. 
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1 Introduction 

The standard model (SM) of particle physics has been extremely successful in describing phe- 
nomena at the highest energies attained thus far. Nevertheless, it is widely believed to be only 
an effective approximation of a more complete theory that supersedes the SM at a higher en- 
ergy scale. Supersymmetry (SUSY) |l]-|8l is generally regarded as one of the likely extensions 
to the SM. The theory is based on the unique way to extend the space-time symmetry group 
underpinning the SM, introducing a relationship between fermions and bosons. 

A low-energy realisation of SUSY, e.g. at the TeV scale, is motivated by the cancellation of the 
quadra tically divergent loop corrections to the Higgs boson mass in the SM (ZUHl. These cor- 
rections are proportional to the masses of the particles that couple to the Higgs boson, where 
the most relevant terms come from the interplay between the masses of the third-generation 
squarks and the largest Yukawa couplings, e.g. of the top quark. In order to avoid a large 
amount of fine-tuning in these loop corrections, the difference in masses between the top quark 
and the third-generation squarks must not be too large 1 9 1 . While the majority of SUSY particles 
(sparticles) may be beyond the reach of the Large Hadron Collider (LHC) at the present beam 
energy and luminosity, the recent discovery of a low-mass Higgs boson candidate |TlOllTT| moti- 
vates "natural" SUSY models in which top and bottom squarks (and gluinos) appear at the TeV 
scale. For R-parity-conserving SUSY [12J, sparticles will be produced in pairs and decay to SM 
particles and the lightest sparticle (LSP), which is generally assumed to be weakly interacting 
and massive. Therefore, the pair production of massive coloured sparticles is expected to result 
in a signature that is rich in jets, especially those originating from bottom quarks, and contains 
a significant amount of missing transverse energy, due to the undetected LSPs. 

This paper summarises an inclusive search for pair production of massive coloured sparticles in 
final states with jets and £j/ performed in pp collisions at a centre-of-mass energy y/s = 8 TeV. 
The analysed data sample corresponds to an integrated luminosity of 11.7 ± 0.5fb _1 |[T3| col- 
lected by the Compact Muon Solenoid (CMS) experiment. Several other searches in this chan- 
nel have been conducted by both the ATLAS and CMS experiments fll4H 25l. The strategy of the 
analysis presented in this paper is based on the kinematic variable ocj, which provides power- 
ful discrimination against multijet production, a manifestation of quantum chromodynamics 
(QCD), while maintaining sensitivity to a wide range of SUSY models. This analysis extends 
previous searches based on a similar strategy with samples of pp collisions at = 7 TeV (23J- 

In order to improve the sensitivity of the analysis to the main production mechanisms of mas- 
sive coloured sparticles at hadron colliders (squark-squark, squark-gluino, and gluino-gluino), 
events with significant £x and two or more energetic jets are categorised according to the jet 
multiplicity. Events with two or three reconstructed jets are used to search for squark-squark 
and squark-gluino production, while events with four or more reconstructed jets probe gluino- 
gluino production. This classification according to the jet multiplicity is a new feature with 
respect to the previous analysis |23|. Moreover, to enhance the sensitivity to third-generation 
squark signatures, events are further categorised according to the number of reconstructed jets 
identified as originating from bottom quarks (b-quark jets). The analysis also considers a large 
dynamic range in the scalar sum of the transverse energies of reconstructed jets in order to 
probe signal models over a large range of mass splittings between the parent sparticle and the 
LSP, including models characterised by a compressed spectrum [26 j. This approach provides 
sensitivity to a wide variety of SUSY event topologies arising from the pair production and de- 
cay of massive coloured sparticles while still maintaining the character of an inclusive search. 
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2 Interpretation with simplified models 

To interpret the results of this search, simplified models ||27H29I are used. These effective mod- 
els include only a limited set of sparticles (production and decay) to enable comprehensive 
studies of individual SUSY event topologies. The result of this search can also be interpreted in 
a range of other relevant models, such as the constrained minimal supersymmetric extension of 
the standard model (CMSSM) |30ti32| or other effective or complete SUSY models that predict 
event topologies with two or more energetic jets and significant J£j- 

In this paper, we focus on the interpretation in two classes of simplified models, the first of 
which describes direct pair-production of squarks, including top and bottom squarks, that de- 
cay to a quark of the same flavour and the LSP. The second class describes gluino-induced 
production of (off-shell) squarks, again including top and bottom squarks, in which gluino- 
gluino pair production is followed by the decay of each gluino to a quark-antiquark pair and 
the LSP. The simplified models considered in this analysis are summarised in Table [T] For each 
model, the LSP is assumed to be the lightest neutralino. 

Table [l]also defines reference models in terms of the parent (gluino or squark) and LSP sparticle 
masses, /n par ent and m^v, respectively, which are used to illustrate potential yields in the signal 
region. In the case of the model D3, a massless LSP is considered. The masses are chosen to be 
reasonably high while still being within the expected sensitivity reach. 

Table 1: A summary of the simplified models considered in this analysis, which involve both 
direct (D) and gluino-induced (G) production of squarks, and their decays. Models Dl and 
G 1 concern the direct or gluino-induced production of first or second-generation squarks only. 
Reference models are also defined in terms of the parent (gluino or squark) and LSP sparticle 
masses. 



Model Production/ decay mode Reference model 

^parent ™LSF 



[GeV] [GeV] 



Dl pp -> qq* -> q^q£? 600 250 

D2 pp -> bb* -> bx^bx® 500 150 

D3 pp -> tt* -4 tj$tx$ 400 

Gl pp -> gg -> qq^qqX? 700 300 

G2 pp -> gg -> bb^bb^ 1 900 500 

G3 pp -> gg -» tt^Q tt^Q 850 250 



3 The CMS apparatus 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam- 
eter, providing a magnetic field of 3.8 T. Within the superconducting solenoid volume are a 
silicon pixel and strip tracker, an electromagnetic calorimeter (ECAL) comprising 75 848 lead- 
tungstate crystals, and a brass/scintillator hadron calorimeter (HCAL). Muons are measured 
in gas-ionisation detectors embedded in the steel flux return yoke of the magnet. Extensive 
forward calorimetry complements the coverage provided by the barrel and endcap detectors. 
The CMS detector is nearly hermetic, which allows for momentum balance measurements in 
the plane transverse to the beam axis. 

CMS uses a right-handed coordinate system, with the origin at the nominal interaction point, 



the x axis pointing to the centre of the LHC ring, the y axis pointing up (perpendicular to the 
plane of the LHC ring), and the z axis along the anticlockwise-beam direction. The polar angle 9 
(radians) is measured from the positive z axis and the azimuthal angle cp (radians) is measured 
in the x-y plane. Pseudorapidity is defined as rj = — ln[tan(0/2)]. 

The silicon pixel and strip tracking systems measure charged particle trajectories with full az- 
imuthal coverage and a pseudorapidity acceptance of 1 77 1 < 2.5. The resolutions on the trans- 
verse momentum (pj) and impact parameter of a charged particle with pj < 40 GeV are typi- 
cally 1% and 15 }im, respectively. Muons are measured in the pseudorapidity range < 2.4. 
Matching muons to tracks measured in the tracking subdetectors results in a pj resolution be- 
tween 1 and 5% for pj < 1 TeV. 

The ECAL has an energy resolution of better than 0.5% for unconverted photons with trans- 
verse energies above 100 GeV. The HCAL, when combined with the ECAL, measures jets with 
a resolution AE/E » 100%/ y/E [GeV] 5%. In the region \rj\ < 1.74, the HCAL cells have 
widths of 0.087 in pseudorapidity and 0.087 in azimuth. In the t]-(p plane, and for \rj\ < 1.48, 
the HCAL cells map onto 5x5 arrays of ECAL crystals to form calorimeter towers project- 
ing radially outwards from close to the nominal interaction point. At larger values of \t]\, the 
size of the towers increases and the matching ECAL arrays contain fewer crystals. Within each 
tower, the energy deposits in ECAL and HCAL cells are summed to define the calorimeter 
tower energies, subsequently used to provide the energies and directions of hadronic jets. 

The first level (LI) of the CMS trigger system, composed of custom hardware processors, uses 
information from the calorimeters and muon detectors to select the most interesting events 
in a fixed time interval of less than 4 ^s. The high-level trigger (HLT) processor farm further 
decreases the event rate, from around 100 kHz to around 300 Hz, before data are stored. 

A more detailed description of the CMS detector can be found in Ref. 1133 1 . 

4 Event reconstruction and selection 
4.1 Definition of ccj 

The ax |25l |34| variable is used to reject multijet events efficiently without significant Ej or 
with transverse energy mismeasurements, while retaining a large sensitivity to new physics 
with final-state signatures containing significant Ej. 

The measurement of Ej typically relies on independent sources of information from each of 
the calorimeter, tracking, and muon subdetectors ||35]| . Relative to other physics objects, this 
measurement is particularly sensitive to the beam conditions and detector performance. This 
difficulty is compounded by the high-energy high-luminosity hadron collider environment at 
the LHC and the lack of precise theoretical predictions for the kinematic properties and cross 
sections of multijet events. 

Given these difficulties, the variable ocj was developed to avoid direct reliance on a measure- 
ment of Ej, instead depending solely on the measurements of the transverse energies and (rel- 
ative) azimuthal angles of jets, which are reconstructed from energy deposits in the calorimeter 
towers II36II . The variable is intrinsically robust against the presence of jet energy mismeasure- 
ments in multijet systems. For dijet events, the otj variable is defined as f2"5ll34"| : 



4 



4 Event reconstruction and selection 



where Ej is the transverse energy of the less energetic jet and Mj is the transverse mass of the 
dijet system, defined as 




where E\, y x x , and p y are, respectively, the transverse energy and x or y components of the 
transverse momentum of jet j ; . 

For a perfectly measured dijet event with ES = E5 and jets back-to-back in (p, and in the limit 
in which the momentum of each jet is large compared with its mass, the value of u.j is 0.5. For 
the case of an imbalance in the measured transverse energies of back-to-back jets, ocj is reduced 
to a value smaller than 0.5, which gives the variable its intrinsic robustness with respect to jet 
energy mismeasurements. Values significantly greater than 0.5 are observed when the two jets 
are not back-to-back and are recoiling against significant, genuine Ej. 

The definition of the ccj variable can be generalised for events with two or more jets as fol- 
lows. The mass scale of the physics processes being probed is characterised by the scalar 
sum of the transverse energy Ej of jets considered in the analysis, defined as Hj = Y%£\ E77 
where n,- e t is the number of jets with Ej above a predefined threshold. The estimator for Ej is 
given by the magnitude of the transverse momenta pj vectorial sum over these jets, defined as 

flj = I fj ]i I . For events with three or more jets, a pseudo-dijet system is formed by com- 
bining the jets in the event into two pseudo-jets. The total Ej for each of the two pseudo-jets 
is calculated as the scalar sum of the measured Ej of the contributing jets. The combination 
chosen is the one that minimises the absolute Ej difference between the two pseudo-jets, AHj. 
This simple clustering criterion provides the best separation between multijet events and events 
with genuine Ej- Equation ^ can therefore be generalised as: 

1 Hj — AHj 1 1 - (AH T /H T ) 

a T = -x ; = - x - v ==L . (3) 

2 ^j Hj 2 _ nj 2 2 Vl - {nj/ H T r 

In the presence of jet energy mismeasurements, the direction and magnitude of the apparent 
missing transverse energy, flj, and energy imbalance of the pseudo-dijet system, AHj, are 
highly correlated. This correlation is much weaker for R-parity-conserving SUSY with each of 
the two decay chains producing the LSP. 

4.2 Physics objects 

Jets are reconstructed from the energy deposits in the calorimeter towers [36J, clustered by the 
infrared-safe ardi-kj algorithm 11371 with a size parameter of 0.5. In this process, the contribu- 
tion from each calorimeter tower is assigned a momentum, the absolute value and the direction 
of which are given by the energy measured in the tower and the position of the tower. The raw 
jet energy is obtained from the sum of the tower energies and the raw jet momentum by the 
vectorial sum of the tower momenta, which results in a nonzero jet mass. The raw jet energies 
are corrected to remove the effects of overlapping pp collisions (pileup) |38ll39 [ and to establish 
a relative uniform response of the calorimeter in t] and a calibrated absolute response in pj. 

The presence of a b-quark jet is inferred from a secondary vertex that is displaced in the trans- 
verse plane with respect to the main interaction vertex. This is achieved with an algorithm |40| 
that incorporates several measurements to build a single discriminating variable that can be 
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used to identify jets originating from bottom quarks with high efficiency and purity. Due to 
the pixel-detector acceptance, b-quark jets are identified only in the region \rj\ < 2.4. In this 
analysis, the discriminator threshold is chosen such that the probability to misidentify (mistag) 
jets originating from light-flavour partons (u, d, s quarks or gluons) as b-quark jets is approxi- 
mately 1% for jets with transverse momenta of 80 GeV [40 1. This threshold results in a b-tagging 
efficiency, i.e. the probability to correctly identify jets as originating from bottom quarks, in the 
range 60-70% EOI , dependent on jet pj. 

The reconstruction of photons, electrons and muons is described below. The presence (or ab- 
sence) of these objects is used to define the event samples for the signal and multiple control 
regions, the latter of which are used to estimate the background contributions from SM pro- 
cesses in the signal region. 

The energy of photons |4"T| is directly obtained from the ECAL measurement, corrected for 
zero-suppression and pileup effects. Various identification criteria must be met in order to cor- 
rectly identify photons with high efficiency and suppress the misidentification of electrons, jets, 
or spurious ECAL noise as photons. These include the requirements that the shower shape of 
the energy deposition in the ECAL be consistent with that expected from a photon, the energy 
detected in the HCAL behind the photon shower must not exceed 5% of the photon energy, and 
no matched hits in the pixel tracker must be found. Isolation from other activity in the event 
is determined through a combination of independent energy sums obtained from each of the 
HCAL, ECAL, and tracker subdetectors within a cone of AR = a/ (A(p) 2 + (At]) 2 = 0.3 around 
the photon trajectory. These sums are corrected for pileup effects and for the contributions from 
the photon itself. 

The energy of electrons [42j is determined from a combination of the track momentum at the 
main interaction vertex, the corresponding ECAL cluster energy, and the energy sum of all 
bremsstrahlung photons attached to the track. Identification criteria similar to those described 
above for photons are applied, with additional requirements on the associated track that con- 
sider the track quality, energy-momentum matching, and compatibility with the main interac- 
tion vertex in terms of the transverse and longitudinal impact parameters. 

The energy of muons |43| is obtained from the corresponding track momentum, combining 
measurements from the muon detectors and both the silicon pixel and strip tracking subdetec- 
tors. Various track quality criteria are considered when identifying muons, as are the transverse 
and longitudinal impact parameters with respect to the main interaction vertex. 

Isolation of muons and electrons is based on a combination of independent sums from the 
HCAL, ECAL, and tracker subdetectors and measured relative to the muon or electron trans- 
verse momentum. The isolation sums are determined for a cone of radius AR = 0.3 (0.4) 
around the electron (muon) trajectory and are corrected for the effects of pileup and for the 
contributions from the lepton itself. 

4.3 Event selection for the signal region 

Events containing non-collision backgrounds are suppressed by requiring at least one vertex of 
high-px tracks to be reconstructed in the luminous region. In the case of multiple vertices, the 
main interaction vertex is defined as the one with the highest scalar sum of pj 2 of all associated 
tracks. 

In order to suppress SM processes with genuine fij from neutrinos in the final state, events are 
vetoed if they contain an isolated electron or muon with pj > 10 GeV. Events with an isolated 
photon with pj > 25 GeV are also vetoed to ensure an all-jet final state. 
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Jets are required to have transverse energy Ej > 50GeV and \tj\ < 3.0. The two highest- 
Ex jets must each satisfy Ej > 100 GeV. These two Ej requirements are relaxed for some 
signal regions, as described below. The highest-Ex jet is additionally required to satisfy \rj\ < 
2.5. Events are vetoed that contain rare, spurious signals from the calorimeters [44J that are 
misidentified as jets. To ensure that the variable Hx is an unbiased estimator of Ej, events are 
vetoed if any additional jet satisfies both Ej > 50 GeV and \t]\ > 3. 

Events are required to have Hj > 275 GeV to ensure high efficiency for the trigger conditions 



used to record the event sample, described in Section 4.4 The signal region is divided into 
eight bins in Hj: two bins of width 50 GeV in the range 275 < Hj < 375 GeV, five bins of 
width 100 GeV in the range 375 < H T < 875 GeV, and a final open bin, Hj > 875 GeV. As in 
Ref . [25], the jet Ej threshold is scaled down to 37 and 43 GeV for the regions 275 < Hj < 325 
and 325 < Hj < 375 GeV, respectively. The threshold for the two highest-Ex jets is also scaled 
accordingly to 73 and 87 GeV. This is done in order to maintain a background composition 
similar to that observed for the higher Hj bins, and also to increase the analysis acceptance for 
SUSY models characterised by compressed spectra. 

Events are further categorised according to the number of jets per event, 2 < nj et < 3 or nj et > 4, 
and the number of reconstructed b-quark jets per event, = 0, 1, 2, 3, or >4. For the category 
of events satisfying nj et > 4 and > 4, the six highest Hj bins are combined to give a final 
open bin of Hj > 375 GeV. 

For events satisfying the selection criteria described above, the multijet background dominates 
over all other SM backgrounds. As discussed in Section |4lj multijet events populate the region 
ax < 0.5. The ax distribution is characterised by a sharp edge at 0.5, beyond which the multijet 
event yield falls by several orders of magnitude. Multijet events with extremely rare but large 
stochastic fluctuations in the calorimetric measurements of jet energies can lead to values of 
olj slightly above 0.5. The edge at 0.5 for multijet events sharpens with increasing Hj, due to 
corresponding increases in the average jet energy and thus an improvement in the jet energy 
resolution. This effect yields an exponential dependence on Hj for the ratio of multijet events 
with a value of ax above and below a given threshold value (larger than 0.5), as described 
further in Section [6l 

The contribution from multijet events is suppressed by many orders of magnitude by requiring 
dj > 0.55. As an example, an event that satisfies both Hj = 275 (875) GeV and ocj = 0.55 must 
also satisfy Hx > 115 (365) GeV. However, certain classes of rare background events can lead 
to values of a.j greater than 0.55, such as those containing beam halo, reconstruction failures, 
spurious detector noise, or event misreconstruction due to detector inefficiencies. These event 
classes, with large, non-physical values of Ej, are rejected by applying dedicated vetoes |35| , 
the most important of which are described below. 

The first example concerns events containing severe energy mismeasurements caused by jet 
trajectories coinciding with inefficient regions in the ECAL (which amount to ~1% of the ECAL 
channel count), or the instrumentation gap between the ECAL barrel and endcap systems at 
= 1.48. These events are identified and vetoed as follows. The negative vector sum of 
jet transverse momenta when jet / is ignored, defined as — Yli=\ i^j Pr\ is determined for each 
ignored jet in turn, 1 < j < Hj et . An azimuthal distance of A(p < 0.5 between the directions 
of jet j and the corresponding vector sum indicates that jet ; has suffered a sufficiently large 
energy mismeasurement to satisfy <xj > 0.55. The event is rejected if the angular distance in the 
(tj, cp) plane between the affected jet and the closest inefficient ECAL region satisfies AR < 0.3. 
Similarly, the event is rejected if the trajectory of the affected jet satisfies Ar/ < 0.3 with respect 
to the ECAL barrel-endcap boundary. 



4.4 Trigger conditions 
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The second example concerns the rare circumstance in which several jets with transverse en- 
ergies below the Ej thresholds and aligned in cp result in significant Ht relative to the value 
of Ej (which is less sensitive to jet Ej thresholds). This type of background, typical of multijet 
events, is suppressed while maintaining high efficiency for SM or new physics processes with 
genuine, significant Ej by requiring flj/Ej < 1.25. The measurement of Ej is provided by the 
particle-flow (PF) reconstruction framework H45ll46l . 

Figure [l] shows the ocj distributions of events with Hj > 375 GeV that satisfy all the selection 
criteria described above except the ocj requirement, categorized according to nj et . An inclusive 
set of trigger conditions is used in order to show the full ocj distribution. The analysis relies 
on data control samples to estimate the contributions from the multijet and non-multijet back- 
grounds, as described in Sections [5] and [6] However, for illustration, the expected yields from 
simulation of multijet events, non-multijet backgrounds with genuine Ej, the sum of these SM 
backgrounds, and an example signal model, are also shown in Fig. [I] The expected yield for 
multijet events that satisfy ocj > 0.55, as given by simulation, is less than ten events and is 
negligible with respect to all other SM backgrounds. Figure [l] highlights the ability of the ccj 
variable to discriminate between multijet events and all other SM or new physics processes 
with genuine Ej in the final state. 
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Figure 1: The ocj distributions of events with Hj > 375 GeV that satisfy all the selection criteria 
described above except the a.j requirement, categorised according to 2 < n^ et < 3 (left) and 
ftjet > 4 (right). An inclusive set of trigger conditions is used to collect the events in data (black 
solid circles with error bars). Expected yields as given by simulation are also shown for multijet 
events (green dash-dotted line), non-multijet backgrounds with genuine Ej as described in 
Section [5] (blue long-dashed line), the sum of all aforementioned SM processes (cyan solid line) 
and the reference signal model D2 (left, red dotted line) or G2 (right, red dotted line). The 
statistical uncertainties for the multijet and SM expectations are represented by the hatched 
areas (visible only for statistically limited bins). The final bin contains all events with ocj > 3. 



4.4 Trigger conditions 

Events are recorded with multiple jet-based trigger conditions, implemented on the HLT com- 
puting farm, that require both Hj and dj to lie above predetermined thresholds, as summarised 
in Table [2] Different trigger conditions are used depending on the analysis Hj bin. The trigger- 
level jet energies are corrected to account for scale and pileup effects. The trigger efficiency, 
defined as the probability with which events that satisfy the signal region selection criteria also 
satisfy the trigger condition, is measured from data for each ftj et category. The efficiency is mea- 
sured using a data sample of }i + jets events recorded by an independent and unbiased trigger 
condition that requires an isolated muon satisfying pj > 24GeV and \t}\ < 2.1. The muon is 
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required to be well separated from the nearest jet by requiring AR(pi,j) > 0.5 and is ignored 
in the calculation of Hj and ccj in order to emulate a genuine JSj signature. 

The thresholds used in the Hj binning scheme are shifted up by 25 GeV with respect to the 
trigger thresholds in order to maintain high efficiency for the Hj component of the trigger 
condition. The measured efficiencies are summarised in Table |2j Non-negligible inefficiencies, 
which are accounted for in the final result, are observed only for the lowest Hj bin. The HLT- 
based trigger conditions are dependent on multiple requirements on quantities determined by 
the LI trigger logic, which include combinations of scalar sums of jet Ej measurements and 
individual Ej thresholds on sub-leading jets. The different efficiencies measured for the two 
nj et categories in the lowest Hj bin are a result of the requirements on LI trigger quantities that 
exhibit non-negligible inefficiencies at very low Hj. 

Table 2: Trigger conditions used to record events for each Hj bin and their efficiencies (with 
statistical uncertainties) measured in data for each Hj bin and n- ]et category. 



Analysis bin Trigger thresholds Trigger efficiency (%) 
Hj [GeV] Hj [GeV] ol t 2 < m et < 3 n jet > 4 



275-325 


250 


0.55 


89 l+°- 4 


03 7+O.6 
-0.6 


325-375 


300 


0.53 


98.7t°J 


98 2 +0 - 4 


375-475 


350 


0.52 


99.0t°j 


QQ 7+0.2 


>475 


400 


0.51 


100.01^ 


lOO.Oj^g 



5 Estimating the non-multijet backgrounds 

5.1 Dominant background processes 

In the absence of a significant contribution from multijet events, the remaining backgrounds in 
the signal region stem from SM processes with significant Ej in the final state. 

For events in which no b-quark jets are identified, the largest backgrounds are from the pro- 
duction of W and Z bosons in association with jets. The decay Z — > vv is the only relevant 
contribution from Z + jets events. For W + jets events, the two relevant sources are leptonic de- 
cays, in which the lepton is not reconstructed or fails the isolation or acceptance requirements, 
and the decay W — > xv in which the T decays hadronically and is identified as a jet. 

For events satisfying > 1, tf production followed by semileptonic decays becomes the most 
important background process. For the subset of events satisfying 1% = 1 and 2 < nj et < 3, 
the total contribution from the W + jets and Z + jets backgrounds is comparable to the tt back- 
ground; otherwise tt production dominates. Events with three or more reconstructed b-quark 
jets originate almost exclusively from tt events, in which one or several jets are misidentified as 
b-quark jets. 

Residual contributions from single-top-quark and diboson production are also expected. 

5.2 Definition of the data control samples 

Three independent data control samples, binned identically to the signal region, are used to 
estimate the contributions from the various background processes. These samples are defined 
by a selection of }i + jets, + jets, and 7 + jets events. The event selection criteria for these 
control samples are defined to ensure that any potential contamination from multijet events 
is negligible. Furthermore, the selections are also expected to suppress contributions from a 
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wide variety of SUSY models (signal contamination) to a negligible level. The selection criteria 
that define the three data control samples are chosen such that the composition of background 
processes and their kinematic properties resemble as closely as possible those of the signal 
region once the muon, dimuon system, or photon are ignored when computing quantities such 
as Hj, AHj, fij, and Oij. This approach emulates the effects, including misreconstruction and 
acceptance, that lead to the presence of these background processes in the signal region. 

The ]i + jets sample is recorded using a trigger condition that requires an isolated muon satis- 
fying pj > 24GeV and \t]\ < 2.1. The event selection requires exactly one isolated muon that 
satisfies stringent quality criteria, pj > 30 GeV, and \rj\ < 2.1 in order for the trigger to be max- 
imally efficient at 88.0 ± 2.0%. Furthermore, the transverse mass of the muon and Ej |45ll46| 
system must be larger than 30 GeV to ensure a sample rich in W bosons. The muon is required 
to be separated from the closest jet in the event by the distance AR > 0.5. The event is rejected if 
two muon candidates are identified that have an invariant mass within a window of ±25 GeV 
around the mass of the Z boson, regardless of the quality and isolation of the second muon 
candidate. No selection requirement on ocj is made in order to increase the statistical precision 
of the predictions derived from this sample, while the impact of relaxing the Uj requirement is 
tested with a dedicated set of closure tests described in Section [5~4l 

The fiji + jets sample uses the same trigger condition and muon identification criteria as the 
p. + jets sample, except that the threshold for the second, lower-px muon is 10 GeV. This leads 
to a trigger efficiency of 95 ± 2%, rising to 98 ± 2% with increasing Hj. The event selection 
requires exactly two oppositely charged, isolated muons satisfying stringent quality criteria, 
and an invariant mass within a window of ±25 GeV around the mass of the Z boson. Both 
muons are required to be separated from their closest jets in the event by the distance AR > 0.5. 
Again, no requirement on ocj is made. 

The 7 + jets sample is selected using a dedicated photon trigger requiring a localised, large 
energy deposit in the ECAL with Ej > 150 GeV that satisfies loose photon identification and 
isolation criteria [41 j. The offline selection requires Hj > 375 GeV, uj > 0.55, and a single 
photon to be reconstructed with Ej > 165 GeV, \r/\ < 1.45, satisfying tight isolation criteria, 
and with a minimum distance to any jet of AR > 1.0. For these selection criteria, the photon 
trigger condition is found to be fully efficient. 

5.3 Method 

The method used to estimate the non-multijet backgrounds in the signal region relies on the 
use of transfer factors, which are constructed per data control sample in bins of Hj, n^ et , and n^. 
These transfer factors are determined from simulated event samples, which are produced as 
follows. The production of W and Z bosons in association with jets is simulated with the M AD- 
GRAPH V5 |47H event generator. The production of tt and single-top quark events is generated 
with POWHEG SHI, and diboson events are produced with PYTHIA 6.4 gU. For all simulated 
samples, PYTHIA 6.4 is used to describe parton showering and hadronisation. The descrip- 
tion of the detector response is implemented using the Geant4 [50J package. The simulated 
samples are normalised using the most accurate cross section calculations currently available, 
usually with next-to-leading-order (NLO) accuracy. To model the effects of pileup, the simu- 
lated events are generated with a nominal distribution of pp interactions per bunch crossing 
and then reweighted to match the pileup distribution as measured in data. 

Each transfer factor is defined as the ratio of expected yields as given by simulation in a given 
bin of the signal region, N^q* 1 , and the corresponding bin of one of the control samples, 
N^q to1 . Each transfer factor is then used to extrapolate from the event yield measured in a 
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data control sample, N^ trol r to an expectation for the event yield in the corresponding bin of 
the signal region, N^f^ 1 , via the expression: 

j. .signal 

^signal _ JV MC y Arcontrol ,a\ 
iV pred - ^control X iV obs ■ W 



Two independent estimates of the irreducible background of Z — > vv + jets events are deter- 
mined from the data control samples comprising Z — > + jets and 7 + jets events, which have 
similar kinematic properties but different acceptances. The Z — > ftp + jets events have similar 
kinematic properties when the two muons are ignored, but a smaller branching fraction, while 
the 7 + jets events have similar kinematic properties when the photon is ignored ||5lJ, but a 
larger production cross section. 

The ji + jets data sample provides an estimate for the total contribution from all other SM 
processes, which is dominated by tt and W-boson production. Residual contributions from 
single-top-quark and diboson production are also estimated. For the category of events satis- 
fying Hj, > 2, in which the contribution from Z — > vv + jets events is suppressed to a negligible 
level, the }i + jets sample is also used to estimate this small contribution rather than using the 
statistically limited + jets and 7 + jets samples. Hence, only the fi + jets sample is used to 
estimate the total SM background for events satisfying n\, > 2, whereas all three data control 
samples are used for events satisfying ti\, < 1. 

In order to maximise sensitivity to potential new physics signatures in final states with multiple 
b-quark jets, a method that improves the statistical power of the predictions from simulation, 
particularly for > 2, is employed. The distribution of is estimated from generator-level 
information contained in the simulation. The number of reconstruction-level jets matched to 
underlying bottom quarks (n^ en ), charm quarks (nf en ), and light-flavoured partons (ft q en ) per 
event, N(n^ en , nf en , n| en ), is recorded in bins of Hj for each ftj et category. The b-tagging effi- 
ciency, e, and mistag probabilities, / c and / q , are also determined from simulation for each Hj 
bin and nj e t category, with each quantity averaged over jet pj and r\. Corrections are applied on 
a jet-by-jet basis to both e, f c , and / q in order to match the corresponding measurements from 
data [40] . This information is sufficient to predict and thus also determine the event yield 
N(ti\,) from simulation for a given Hj bin and ftj et category with the expression: 

N(n h ) = E E ( N(nf n , nf n , «f n ) x P b x P c x P q ) , (5) 

"jet n b 

where nj^ 8 , n' ag , and n q g are the number of times that a reconstructed b-quark jet is identified 
as originating from an underlying bottom quark, charm quark, or light-flavoured parton, re- 
spectively, and Pb = P( n b S ' n b en ' £ )' = P( n c aS ', n c en ,fc), and P q = P(n q s ;n q en ,/ q ) are the 
binomial probabilities for this to happen. The outer summation considers all possible combi- 
nations of w gen , nf en , and n q en that satisfy ftj et = n gen + nf en + n q en , while the inner summation 

considers all possible combinations of n£ s , n' a8 , and n q 8 that satisfy n b = + n' as + n q s . 

The predicted yields are found to be in good statistical agreement with the yields obtained 
directly from the simulation in the bins with a significant population. The method exploits 
the ability to make precise measurements of N(n gen , nf en ,n q en ), e, f c , and / q independently of 
n\„ which means that event yields for a given b-quark jet multiplicity can be predicted with a 
higher statistical precision than obtained directly from simulation. Precise measurements of f c 
and fa are particularly important for events with > 3, which often occur in the SM because 
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of the presence of mistagged jets in the event. In this case, the largest background is tt, with two 
correctly tagged b-quark jets and an additional mistagged jet originating from a charm quark 
or light-flavoured parton. 

5.4 Systematic uncertainties on transfer factors 



As described in Section 5.3 the method to estimate the background contributions from SM 
processes with significant £t is based on an extrapolation from a measurement in a control 
sample to a yield expectation in the signal region. This approach aims to minimise the sensi- 
tivity to simulation mismodelling, as many systematic biases are expected largely to cancel in 
the ratios used to define the transfer factors. However, a systematic uncertainty is assigned to 
each transfer factor to account for theoretical uncertainties [|5T] and residual biases in the sim- 
ulation modelling of kinematics (e.g. acceptances) and instrumental effects (e.g. reconstruction 
inefficiencies). 

The magnitudes of the systematic uncertainties assigned to the transfer factors are determined 
from a representative set of closure tests in data. These tests use yields from an event category 
in one of the three independent data control samples, along with the corresponding transfer 
factors obtained from simulation, to predict the yields in another event category or data control 
sample following the prescription defined in Eq. Q. As stated previously, the contamination 
from multijet events or any potential signal is expected to be negligible. Therefore, the closure 
tests carried out between control samples probe the properties of the relevant SM non-multijet 
backgrounds. 

Thirteen sets of closure tests are chosen to probe key ingredients of the simulation modelling 
that may introduce biases in the transfer factors. Each set comprises up to eight independent 
tests in bins of Hj. Five sets of closure tests are performed independently for each of the two 
ftjet categories, and a further three sets are common to both categories, as shown in Fig. [2] For 
each nj e t category, the first three sets of closure tests are carried out within the ]i + jets sample, 
and probe the modelling of the aj distribution in genuine f.j events (circles), the relative com- 
position between W + jets and top events (squares), and the modelling of the reconstruction 
of b-quark jets (triangles), respectively. The fourth set (crosses) addresses the modelling of the 
vector boson samples by connecting the }i + jets and + jets control samples, with the former 
sample rich in W + jets events (and also with a significant contribution from top events) and 
the latter in Z + jets events. The fifth set (solid bullets) deals with the consistency between 
the Z — > }i}i + jets and 7 + jets samples, which are both used to provide an estimate of the 
Z — > vv + jets background. Three further sets of closure tests (inverted triangles, diamonds, 
asterisks), one per data control sample, probe the simulation modelling of the n- ]et distribution. 

All sets of closure tests demonstrate, given the statistical precision of each test, that there are 
no significant biases or dependencies on Hj exhibited by the transfer factors obtained from 
simulation. Table [3] summarises the results obtained from constant and linear polynomial fits 
to each set of closure tests for the two nj et categories. The table also lists the best fit values 
and uncertainties for the constant polynomial fits, which indicate the level of closure averaged 
across the full Hj range considered in the analysis. All tests are either statistically compatible 
with zero bias or at the level of a few percent or less. Finally, Table [3] also summarises the 
best fit values of the slopes of the linear polynomial fits, which are typically of the order 10~ 4 , 
corresponding to a percent-level change per 100 GeV. However, in all cases, the best fit values 
are fully compatible with zero, indicating that the level of closure is Hj-independent. The x 2 
and number of degrees of freedom (d.o.f.) of each fit are also quoted and indicate a reasonable 
goodness-of-fit in all cases except one, which concerns the simulation modelling of the Hj e t 
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Figure 2: Sets of closure tests that probe for possible Hj-dependent biases associated with the 
transfer factors obtained from simulation, for the two event categories satisfying 2 < nj et < 3 
(left) and n- ]et > 4 (right). Also shown are shaded bands that represent Hj-dependent system- 
atic uncertainties. 



distribution in the }i + jets sample. The large x 2 value is mainly attributable to a single outlier 
in the bin 675 < Hj < 775 GeV rather than any significant trend in Hj. 

Table 3: Results from constant and linear polynomial fits to sets of closure tests performed for 
each ftj e t category. The symbol identifies the set of closure tests in Fig. [2] The final four rows 
probe the simulation modelling of the Hj et distribution. The + indicates the fit repeated with a 
single outlier removed. 

Constant polynomial fit Linear polynomial fit 



"jet 


Symbol 


Set of closure tests 


Constant 


X 2 /d.o.f. 


Slope 
(10- 4 GeV" 1 ) 


X 2 /d.o.f. 


2-3 


o 


a T < 0.55 -> a T > 0.55 (ff + jets) 


-0.06 ± 0.02 


2.43/7 


-1.3 ±2.2 


2.10/6 


2-3 


□ 


b tags — > 1 b tag(f( + jets) 


0.07 ± 0.02 


1.49/7 


-1.6 ±1.6 


0.54/6 


2-3 


A 


1 b tag — > 2 b tags(f( + jets) 


-0.07 ±0.03 


4.19/7 


-2.7 ±3.0 


3.41/6 


2-3 


X 


}i + jets — > }t}i + jets 


0.10 ±0.03 


5.64/7 


-1.1 ±2.3 


5.40/6 


2-3 


• 


fifi +jets — >■ 7 + jets 


-0.06 ± 0.04 


5.93/5 


4.2 ± 4.3 


4.98/4 


>4 





txj < 0.55 -> a T > 0.55 (f( + jets) 


-0.05 ± 0.03 


9.58/7 


3.0 ± 2.9 


8.47/6 


>4 


□ 


b tags — > 1 b tag(f( + jets) 


-0.03 ± 0.03 


5.88/7 


-1.0 ±1.9 


5.59/6 


>4 


A 


1 b tag — ► 2 b tags(f( + jets) 


-0.02 ± 0.03 


7.35/7 


1.1 ±2.2 


7.08/6 


>4 


X 


}i + jets — )• jiji + jets 


0.08 ± 0.07 


12.9/7 


4.8 ± 4.3 


11.7/6 


>4 


• 


+jets — >■ 7 + jets 


-0.03 ±0.10 


2.85/5 


-4.0 ±7.0 


2.52/4 


>2 


V 


2 < «jet < 3 -> Mjet > 4(fi + jets) 


-0.03 ± 0.02 


17.3/7 


0.0 ±1.0 


17.3/6 


>2 + 


V 


2 < «jet < 3 -S- Wjet > 4(ft + jets) 


-0.04 ±0.01 


6.10/6 


-1.4 ±1.1 


4.46/5 


>2 





2 < Hjet < 3 -> wjet > 4(7 + jets) 


0.12 ±0.05 


2.42/5 


6.0 ±4.7 


0.77/4 


>2 


* 


2 < Mj et < 3 — > Mjet > 4(f(^ +jets) 


-0.04 ±0.07 


9.76/7 


4.9 ± 4.4 


8.51/6 



Once it is established that no significant bias or trend is observed for any set of closure tests, un- 
corrected systematic uncertainties on the transfer factors are determined for five independent 
regions in Hj: 275-325, 325-375, 375-575, 575-775, and > 775 GeV. Conservative estimates for 
the systematic uncertainties are based on the variance in the level of closure for all individ- 
ual tests, weighted according to the statistical uncertainties associated with each test, within 
a given H T region. This procedure yields estimates of 10% (10%), 10% (10%), 10% (10%), 20% 
(20%), and 20% (30%) for the five Hj regions defined above for events satisfying 2 < n; et < 3 
(ftjet > 4), as indicated by the shaded bands in Fig. [2] 

The effect on the transfer factors of uncertainties related to the modelling of b-quark jets in 
simulation is found to be negligible, at the sub-percent level, in comparison to the aforemen- 
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tioned Hj-dependent systematic uncertainties. Hence, the same estimates for the systematic 
uncertainties are used for each category which are treated as uncorrelated between the n\, 
categories. 



6 Estimating the multijet background 

The contribution from multijet events is expected to be negligible, at or below the percent- 
level relative to the yields expected from non-multijet backgrounds, even for the most inclusive 
definition of the signal region, defined by Hj > 275 GeV, ax > 0.55, and no requirement on 
njet or Hi,. The expected yield is further suppressed to <C 1 event with the application of more 
stringent thresholds on any of the variables Hj, n- ]et , or n^. 

Any potential contamination from multijet events can be estimated by exploiting the Hj de- 
pendence of the ratio of events with a value of ocj above and below some threshold, R Kj (Hj). 
This dependence on Hj is modelled as a falling exponential function, R Kj (Hj) = Ae~ kHj p5| , 
where the parameters A and k are the normalisation and decay constant parameters, respec- 
tively. The exponential model is validated in a multijet-enriched data sideband, defined by the 



event selection criteria for the signal region, described in Section 4.3 but with the requirement 
Ht/^t > 1.25. A measurement of the decay constant k is made in a further multijet-enriched 
sample defined by the event selection criteria for the signal region but with the requirement 
olj < 0.55. 

The estimate of the multijet contamination in the signal region for a given Hj bin is determined 
from the product of the ratio R aT , as given by the exponential model, and the yield in a data 
control sample defined by the event selection for the signal region but with the requirement 
dlj < 0.55. This event sample is recorded with a set of trigger conditions that require only Hj 
to be above the same thresholds as used by the signal region triggers listed in Table [2] 

Further details on the exponential model and its use in the likelihood model are found in Sec- 
tion 

7 Confronting data with the SM-only hypothesis 

For a given category of events satisfying requirements on both ftj et and n^, a likelihood model 
of the observations in multiple data samples is used to obtain a consistent prediction of the SM 
backgrounds and to test for the presence of a variety of signal models. It is written as: 



"y r 



(0 < n h < 1) 
(n b > 2) 



(6) 
(7) 



where L^r describes the yields in the eight Hj bins of the signal region where exactly b- 
quark jets are required. In each bin of Hj, the observation is modelled as a Poisson-distributed 
variable about the sum of the SM expectation and a potential signal contribution (assumed to 
be zero in the following discussion), where the SM expectation is the sum of the multijet and 
non-multijet components. The non-multijet component is related to the expected yields in the 
\i + jets, + jets, and 7 + jets control samples via the transfer factors derived from simulation, 

The likelihood functions L„, L }lfl , and L 7 describe the yields in 



as described in Section 5.3 



the Hj bins of the pi + jets, + jets, and 7 + jets control samples. For the category of events 
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satisfying n\, > 2, only the ji + jets control sample is used in the likelihood to determine the total 
contribution from all non-multijet SM backgrounds in the signal region. The estimate of the 
contribution from multijet events in a given Hj bin of the signal region relies on the exponential 
model R Kj (Hj) = Ae~ kHj , as described in Section^ The systematic uncertainties (10-30%) 
associated with the SM background estimates are accounted for with nuisance parameters in 
the likelihood function, the measurements of which are treated are assumed to follow a log- 
normal distribution. 

In order to test the compatibility of the observed yields with the expectations from only SM 
processes, the likelihood function is maximised over all fit parameters. Table |4] summarises 
the observed yields and the SM expectations in bins of Hj for events in the signal region cate- 
gorised according to ftj et and n\,. No significant excess above the SM expectation is observed in 
the signal region and the control samples are well described by the SM hypothesis. For each of 
the eight categories of events defined by requirements on m e t and nj-,, the goodness-of-fit of the 
SM-only hypothesis is determined by considering simultaneously up to eight bins in Hj from 
the signal region and up to 22 bins from the three control samples. The p-values obtained are 
found to be uniformly distributed in the range 0.0-1.0. 

Comparisons of the observed yields and the SM expectations in bins of Hj for events cate- 
gorised according to Hj et and containing exactly zero, one, or two b-quark jets are shown in 
Figs. |3j|4j and [5} respectively. Similarly, Fig. [6] shows the Hj-binned observed yields and SM 
expectations for events satisfying n^ et > 4 and = 3 (left) or > 4 (right). For illustra- 
tion, Figs.[3]-[6]include the expected yields from various reference models, as defined in Table]!] 
Figure [7] (left column) shows the observed yields and SM expectations in the Hj bins of the 
}i + jets, + jets, and 7 + jets control samples for events satisfying 2 < nj e t < 3 and = 0. 
Figure [7] (right column) shows the observed yields and SM expectations in the Hj bins of the 
}i + jets sample for events satisfying n ]et > 4 and n\, = 2,n\ y = 3, or > 4. 

The maximum-likelihood values for the decay constant and normalisation parameters, k and 
A, of the exponential model for the multijet background are obtained independently for each of 
the eight event categories. The value of the nuisance parameter k is constrained via a measure- 
ment in a multijet-enriched data sideband, as described in Section|6] No constraint is applied to 
the normalisation term. In the nominal fit, the maximum-likelihood value of the normalisation 
parameter for each event category is found to be compatible with zero within uncertainties. 
Furthermore, the expected yields obtained from an alternate fit, in which the normalisation 
parameters are fixed to zero, agree well with those obtained from the nominal fit. 

8 Interpretations of the results 

Limits are set in the parent sparticle and LSP mass parameter space of the simplified models 
listed in Table [l] The CLs method |52l |53l is used to compute the limits, with the one-sided 
(LHC-style) profile likelihood ratio as the test statistic |54]. The sampling distributions for 
the test statistic are built by generating pseudo-data from the likelihood function, using the 
respective maximum-likelihood values of the nuisance parameters under the SM background- 
only and signal-plus-background hypotheses. Signal contributions in each of the data samples 
are considered, though the only significant contribution occurs in the signal region and not the 
control samples. Table [5] specifies the event categories, defined in terms of nj et and n^, used to 
provide interpretations in the different simplified models. 

Event samples for the simplified models are generated at leading order with pythia 6.4 [49 1. 
Inclusive, process-dependent, NLO calculations of SUSY production cross sections, with next- 
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Table 4: Event yields observed in data and SM expectations with their associated uncertainties 
in bins of Hj for events in the signal region that are categorised according to ftj et and n\,. The 
final Hj > 375 GeV bin is inclusive for the n- jet > 4 and n\, > 4 category. 
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Figure 3: Event yields observed in data (solid circles) and SM expectations with their associated 
uncertainties (solid lines with bands) in bins of Hj for the signal region when requiring exactly 
zero b-quark jets and 2 < n- jet < 3 (left) or Mj et > 4 (right). For illustration only the expectations 
for the reference mass points of the signal models Dl (left, red dashed line) and Gl (right, red 
dashed line) are superimposed on the SM-only expectations. 



Table 5: A summary of the event categories used to provide an interpretation in the various 
simplified models considered in this paper. 
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Figure 4: As for Fig. |3j but requiring exactly one b-quark jet and 2 < n^ et < 3 (left) or nj et > 4 
(right). Example signal yields are for the reference mass points of the signal models D2 (left, 
red dashed line) and D3 (right, red dashed line). 
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Figure 5: As for Fig. |3j but requiring exactly two b-quark jets and 2 < n ]et < 3 (left) or Wj et > 4 
(right). Example signal yields are for the reference mass points of the signal models D2 (left, 
red dashed line) and D3 (right, red dashed line). 



'cms, 11.7 fb" 1 , is = 8 TeV 
— I — Data (signal region, n b = 3; n > 4) 

Standard Model ± Expected Unc. 

SM + pp^gg,g->bb"jc° 

(m_= 900 GeV, m „ = 500 GeV) 




300 400 500 600 700 



800 



c 10' 



10 



900 
Ht (GeV) 



10" 



CMS, 11.7 fb" 1 , is = 8 TeV 

— + — Data (signal region, n b > 4; n. > 4) 

Standard Model ± Expected Unc. 

SM + pp->gg,g->tfx° 

(m.= 850 GeV, m = 250 GeV) 



280 300 320 340 360 380 400 420 440 460 

hir (GeV) 



Figure 6: As for Fig. |3j but requiring nj et > 4 and exactly three (left) or at least four (right) 
b-quark jets. Example signal yields are for the reference mass points of the signal models G2 
(left, red dashed line) and G3 (right, red dashed line). 
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Figure 7: Event yields observed in data (solid circles) and SM expectations with their associated 
uncertainties (solid lines with bands) in bins of Hj for: the ]i + jets (top left), + jets (middle 
left), and 7 + jets (bottom left) control samples when requiring 2 < nj et < 3 and exactly zero b- 
quark jets; and the }i + jets control sample when requiring nj et > 4 and exactly two (top right), 
three (middle right), or at least four (bottom right) b-quark jets. 
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8 Interpretations of the results 



to-leading-logarithmic (NLL) corrections, are obtained with the program PROSPINO |55l - l60| . 
The samples are generated using the CTEQ6L1 [61 j PDFs. The distribution of the number of 
pp interactions per bunch crossing for the simulated samples matches that observed in data. 

Various experimental uncertainties on the expected signal yield are considered for each inter- 
pretation. Signal acceptance in the kinematic region defined by < m pa rent — wlsp < 175 GeV 
or wzp aren t < 300 GeV is due in part to the presence of initial-state radiation. Given the large as- 
sociated uncertainties from simulation for this kinematic region, no interpretation is provided. 
Otherwise, the experimental systematic uncertainties are determined for each point in the mass 
parameter space of each simplified model. Models are categorised according to the mass split- 
ting between the parent sparticle and the LSP, with those satisfying 175 < m paient — Wlsp < 
350 GeV deemed to be characterised by a compressed spectrum. For a given category of model, 
i.e. with a compressed spectrum or otherwise (as defined above), the point-by-point variations 
in the mass parameter space are small with respect to the absolute values of the systematic 
uncertainties. 

Estimates of the various systematic uncertainties for models with a compressed spectrum or 
otherwise are summarised in Tables [6] and |7j respectively. Contributions from the analysis se- 
lection are dominated by uncertainties on the PDFs, jet energy scale (JES), and modelling of 
the efficiency and mistag probability of b-quark jets in simulation. The total systematic un- 
certainties provided in the tables also account for the uncertainty of 4.4% on the luminosity 
measurement [ 13 1 and contributions from other event selection criteria, such as: the trigger 
conditions; the removal of events containing isolated muons, electrons, or photons; and filters 



to suppress classes of rare, pathological events, as described in Section 4.3 Each of these in- 
dividual contributions is below 4%. The total systematic uncertainty on the expected signal 
yield for the various simplified models is found to be in the range 12%-23% and is accounted 
for with a nuisance parameter, the measurement of which is assumed to follow a lognormal 
distribution. 

Table 6: Estimates of the dominant systematic uncertainties (%), defined in the text, on the 
analysis efficiency for various simplified models that are characterised by a small mass split- 
ting (i.e. compressed spectrum) between the parent sparticle and LSP. The totals also re- 
flect contributions from additional systematic uncertainties described in the text. The region 
^parent — w?lsp < 350 GeV is kinematically forbidden for the G3 model. 



Model 


Dl 


D2 


D3 


Gl 


G2 G3 


PDF 


10.0 


10.0 


10.0 


10.0 


10.0 


JES 


4.1 


4.8 


6.5 


5.6 


7.3 


b-tagging 


2.4 


2.2 


0.8 


3.1 


2.7 


Total 


12.9 


13.1 


13.9 


13.9 


14.5 



Table 7: Estimates of the dominant systematic uncertainties (%), defined in the text, on the 
analysis efficiency for various simplified models that are characterised by a large mass split- 
ting between the parent sparticle and LSP. The totals also reflect contributions from additional 
systematic uncertainties described in the text. 



Model 


Dl 


D2 


D3 


Gl 


G2 


G3 


PDF 


10.0 


10.0 


10.0 


10.0 


10.0 


10.0 


JES 


1.1 


0.9 


3.5 


0.8 


1.5 


0.5 


b-tagging 


5.8 


2.7 


1.6 


6.6 


10.1 


19.4 


Total 


13.4 


12.3 


12.9 


14.0 


16.0 


23.0 



Figure [8] shows the observed upper limit on the production cross section at 95% confidence 
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level (CL) as a function of the parent sparticle and LSP masses for various simplified models. 
The point-to-point fluctuations are due to the finite number of pseudo-experiments used to 
determine the observed upper limit. Also indicated are the observed excluded region assum- 
ing NLO+NLL cross sections for squark pair production in the limit of very massive gluinos 
(or vice versa), the observed excluded region when varying the production cross section by 
its theoretical uncertainty and the expected excluded region with the ±1 standard-deviation 
variations. 

Table [8] lists the expected signal yields and analysis efficiencies in the region Hj > 375 GeV for 
each of the reference models defined in Table]!] The yields and efficiencies are summed over the 
individual event categories used for each interpretation, as listed in Table [5] The observed and 
expected upper limits (95% CL) on the cross section are also quoted, which can be compared 
with the NLO+NLL SUSY production cross section and its theoretical uncertainty. 

Table 8: Summary of expected yields, analysis efficiencies, and upper limits for the reference 
models defined in Tablejl] using the event categories defined in Table [5] The first row specifies 
the reference model. The second and third rows quote the expected yield and analysis effi- 
ciency (with statistical uncertainties) for the region Hj > 375 GeV. The fourth row quotes the 
NLO+NLL SUSY production cross section (with theoretical uncertainty), and the fifth and sixth 
rows quote the observed and expected upper limits (95% CL) on the production cross section. 



Reference model 


Dl 


D2 


D3 


Gl 


G2 


G3 


Expected yield 


358.3 + 8.9 


78.1 + 2.4 


90.6 + 2.4 


416 ± 13 


52.0+1.7 


25.3 ± 0.7 


Efficiency (%) 


16.0 + 0.4 


10.2 + 0.3 


2.9 + 0.1 


10.4 + 0.3 


9.4 + 0.3 


2.9 + 0.1 


Theoretical cross section [fb] 


245 ± 44 


86 + 13 


357 + 51 


434 ± 81 


60+14 


97 + 21 


Observed upper limit [fb] 


113.2 


42.3 


360.8 


103.0 


15.0 


46.2 


Expected upper limit [fb] 


103.1 


31.2 


240.6 


65.2 


12.3 


35.3 



The estimates of mass limits are determined from the observed exclusion based on the theo- 
retical production cross section, less one-standard-deviation uncertainty. The most stringent 
mass limit on the parent sparticle, ?«p|®g nt , is generally obtained at low LSP masses. This limit 
weakens for larger values of LSP mass until a value »z£|p * s reacne d/ beyond which no limit can 
be set. Generally speaking, the excluded mass range for m pa rent is bounded from below by the 
kinematic region considered for each model, up to Wlsp ~ m LSP- This implies that the value 
of mk|p is relatively insensitive to the mass splitting, and is generally valid for models charac- 
terised by compressed spectra or otherwise. Whether this statement holds true for very small 
mass splittings, satisfying m patent — m^sp < 175 GeV, requires further studies of systematic un- 
certainties on signal acceptance such as those due to, for example, JES or the identification of 
b-quark jets. These features can be illustrated in Fig. |8j (middle), where ^parent — 

950 GeV for 

the region wjlsp < 350 GeV. This limit weakens to 900 GeV when the LSP mass reaches 425 GeV. 
However, for LSP masses above 450 GeV, no mass range can be excluded for gluinos decaying 
to first- or second-generation quarks. The reach in LSP mass is stable down to mass splittings 
Of Wparent - m LSP « 175 GeV. 

Table ^summarises the most stringent observed and expected mass limits, in terms of m^f^ 

and wz£gp, obtained for the simplified models considered in this paper. The observed exclusion 
for each simplified model is generally weaker than expected at the level of 1-2 standard devi- 
ations. This feature is attributed to the small upward fluctuations in data in either the region 
Hj > 875 GeV for the = category or 475 < Hj < 675 GeV for the categories of events sat- 
isfying 1 < 1% < 2. Candidate events in these regions have been examined and do not exhibit 
any non-physical behaviour. The expected search sensitivity has improved with respect to the 
analysis based on the \/s = 7TeV dataset |23| by as much as 225 and 150 GeV for and 
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Figure 8: Observed upper limit on the production cross section at 95% CL (indicated by the 
colour scale) as a function of the parent and LSP sparticle masses for simplified models involv- 
ing: the direct production of first- and second-generation squarks (Dl, top left) and bottom 
squarks (D2, top right); and pair-produced gluinos each decaying to the LSP and pairs of light 
quarks (Gl, middle), bottom quarks (G2, bottom left), or top quarks (G3, bottom right). The 
black solid thick line indicates the observed exclusion assuming NLO+NLL SUSY production 
cross section. The black solid thin lines represent the observed exclusions when varying the 
cross section by its theoretical uncertainty. The purple dashed thick (thin) line indicates the 
median (±lc") expected exclusion. 
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m^lp, respectively. 

Table 9: Summary of the mass limits obtained for various simplified models. The limits indicate 
the observed (expected) search sensitivity for each simplified model, where m^ ent and m^|p 
represent the largest mass beyond which no limit can be set for squarks or gluinos and the LSP, 
respectively. No observed exclusion is possible in the mass parameter space considered for the 
model D3. 



Model 


Dl 


D2 


D3 


Gl 


G2 


G3 


<£nt [GeV] 
m^fpt [GeV] 


775 (850) 
325 (350) 


600 (675) 
200 (250) 


- (520) 
- (100) 


950 (1050) 
450 (550) 


1125 (1200) 
650 (700) 


950 (1075) 
325 (375) 



Figure [9] shows the observed upper limit at 95% CL on the production cross section as a func- 
tion of the top-squark mass (m -j) for the model D 3 when considering different LSP masses in the 
range 0-150 GeV. No exclusion on possible top-squark masses is observed when considering 
the theoretical production cross section, less lc uncertainty. However, the expected exclusion 
covers the ranges 300-520, 320-520, and 420-480 GeV for m LS p = 0, 50, and 100 GeV, respec- 
tively. No exclusion is expected for the LSP with a mass greater than 100 GeV. The expected 
reach for the D3 model is summarised in Table |U 



9 Summary 

An inclusive search for supersymmetry with the CMS experiment is reported, based on a data 
sample of pp collisions collected at y/s = 8 TeV, corresponding to an integrated luminosity of 
11.7 ± 0.5 fb _1 . Final states with two or more energetic jets and significant JS-r, as expected from 
the production and decay of massive squarks and gluinos, have been analysed. 

The analysis strategy is to maximise the sensitivity of the search to a wide variety of SUSY 
event topologies arising from squark-squark, squark-gluino, and gluino-gluino production and 
decay, particularly those with third-generation squark signatures, while still maintaining the 
inclusive nature of the search. The signal region is binned according to the number of recon- 
structed jets, the scalar sum of the transverse energy of jets, and the number of jets identified 
to originate from bottom quarks. The sum of standard model backgrounds per bin has been 
estimated from a simultaneous binned likelihood fit to event yields in the signal region and 
}i + jets, fiji + jets, and 7 + jets control samples. The observed yields in the signal region are 
found to be in agreement with the expected contributions from standard model processes. Lim- 
its are set in the SUSY particle mass parameter space of simplified models, with an emphasis 
on the different production mechanisms of coloured SUSY particles, third-generation squark 
signatures, and compressed-spectrum scenarios. The results can also be used to perform inter- 
pretations in other relevant models, such as the CMSSM. 

In the context of simplified models, gluino masses below ~1 TeV are excluded at the 95% CL 
under the assumptions that gluinos are produced in pairs and each decays to a quark-antiquark 
pair and a light LSP via an off-shell squark. The mass limit varies in the range 950-1125 GeV 
depending on the squark flavour. The most constraining mass limits on the gluino-induced LSP 
are in the range 325-650 GeV depending on the decay mode. For the direct pair production of 
squarks, each of which is assumed to decay to a quark and light LSP, masses below 775 and 
600 GeV are excluded (95% CL) for first- and second-generation squarks and bottom squarks, 
respectively. The most constraining mass limits on the LSP from the decays of first- and second- 
generation squarks and bottom squarks are 325 and 200 GeV, respectively. No exclusion is 
observed for the direct pair production of top squarks, each of which is assumed to decay to 
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Figure 9: Excluded cross sections versus top-squark mass for the model D3, in which pair- 
produced top squarks each decay to a top quark and the LSP with a mass wJlsp = (top left), 50 
(top right), 100 (bottom left), and 150 GeV (bottom right). The observed upper limit (95% CL) 
on the production cross section is shown as a function of m~ t (solid line), along with the expected 
upper limit and ±lc experimental uncertainties (long-dashed line with shaded band), and the 
NLO+NLL top-squark pair-production cross section and theoretical uncertainties (dotted line 
with shaded band). 
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a top quark and the LSP. However, an exclusion is expected for top squark masses as high as 
~500 GeV and an LSP mass as high as 100 GeV. The limits on the LSP masses are also generally 
valid for compressed-spectrum models with mass splittings between the parent sparticle and 
LSP as low as ~200 GeV. 

The analysis strategy reported here, in conjunction with the increase in centre-of-mass energy 
to 8 TeV, has increased the coverage of SUSY parameter space with respect to previous searches. 
However, a large range of SUSY parameter space still remains to be probed by the LHC. 
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